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Standard Monte Carlo simulations have been performed for 3-dimensional U(l) lattice gauge 
model on improved lattices to measure the wave functions and the size of the scalar and tensor 
glueballs. Our results show the radii of ~ 0.60 and ~ 1.12 in the units of string tension, or ~ 0.28 
and ~ 0.52 fm, for the scalar and tensor glueballs, respectively. At finite temperature we see 
clear evidence of the deconfined phase, and the transition appears to be similar to that of the 
two-dimensional XY model as expected from universality arguments. Preliminary results show no 
significant changes in the glueball wave functions and masses in the deconfined phase. 



PACS numbers: 11.15.Ha, 12.38.Gc,11.15.Me 



I. INTRODUCTION 



The prediction of glueball masses has longbeen at- 
tempted in lattice gauge theory calculations 0J2>EIE!>O1- 
These calculations show that the lowest-lying scalar, ten- 
sor and axial vector glueballs lie in the mass region of 
1-2.5 GeV. While there is a long history of glueball mass 
calculations in lattice QCD, little is known about the 
glueballs besides their masses . Accurate lattice calcu- 
lations of their size, matrix elements and form factors 
would help considerably in their experimental identifica- 
tion. 

Glueball wave functions and sizes have been studied 
m the past H 0, H 0, but much of the early work 
contains uncontrolled systematic errors, most notably 
from discretisation effects. The scalar glueball is par- 
ticularly susceptible to such errors for the Wilson gauge 
action, due to the presence of a critical end point of a 
line of phase transitions (not corresponding to any phys- 
ical transition found in QCD) in the fundamental-adjoint 
coupling plane. As this critical end-point (which defines 
the continuum limit of a </> 4 scalar field theory) is neared, 
the coherence length in the scalar channel becomes large, 
which means that the mass gap in this channel becomes 
small; glueballs in other channels seem to be affected very 
little. Results in which the scalar glueball was found to 
be significantly smaller than the tensor were most likely 
due to contamination of the scalar glueball from this non- 
QCD critical point pfl. On the other hand, the calcula- 
tions using operator overlaps obtained from variational 
optimization for improved lattice gauge action, which are 
designed to avoid spurious critical point, show that the 
scalar and tensor glueballs were of typical hadronic di- 
mensions pi . A straightforward procedure to address 
the controversy over glueball size is to measure the glue- 



ball wave function, much in the same way as the meson 
and baryon wave functions were measured jlCj . 

As a first test of our method we study the low-lying 
scalar and tensor glueballs and their wave functions for 
the 3-dimensional U(l) model with tadpole improved 
Symanzik gauge action. This theory is perhaps the best 
understood of all lattice gauge theories at zero tempera- 
ture, which makes it testing ground for new analytic and 
numerical methods. The relevance of the model to QCD 
at finite temperature 111), including confinement, chiral 
symmetry breaking |l2j and existence of a mass gap^l, 
suggest that a comparison of the respective mass spec- 
tra should be informative. Our techniques for calculating 
the glueball wave functions from Wilson loop operators 
are outlined in Sec. [HI We present and discuss our re- 
sults at zero temperature in Sec. IIIII We extend our 
method to examine the wave functions and masses at fi- 
nite temperature in this section. Here we briefly discuss 
an interpretation of the deconfincmcnt in terms of the 
power-law behaviour of the correlation function. Section 
IIVI is devoted to the summary and concluding remarks. 



II. WAVE FUNCTIONS OF GLUEBALLS 

In contrast with the techniques used in Ref. 0, we 
measure our lattice operators from spatially connected 
Wilson loops. Glueballs are colour singlet states and one 
should be able to construct them with closed-loop paths 
which are gauge invariant. The choice of such loops elim- 
inate the need for gauge fixing 1 . Although the calcula- 
tions in Ref. have produced some interesting results, 
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1 It should be noted that gauge-invariant Wilson loops have a 4 de- 
pendence compared to a 2 dependence of two-link operator used 
in Ref. []]. The lower dimension operator yields a linear depen- 
dence in the correlation function as opposed to a 5 dependence for 
Wilson loops. This improves the glueball signal as a is reduced. 
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the approach suffers from a basic problem: the observ- 
ables are calculated from a lattice version of the 2-glue 
operator, which risks a mixture of glueball states with 
flux states 2 . 

In this paper we take a more direct approach to the 
problem. We measure the observables in a three-step 
procedure. First, we calculate the lattice operator 



$(f,f)=^[#M) 



(1) 



where <f> is the plaquette operator and $ measures the 
two plaquette or two-loop component of the glueball wave 
function. The r dependence will be reflected in the length 
of links required to close the loops. From a suitable lin- 
ear combinations of rotation, parity inversions and real 
or imaginary parts of the operators involved in <E>, one 
can construct glueball operators with desired quantum 
numbers. Since we want to explore the nature of wave 
functions, we focus only on the low- lying "symmetric" 
and "antisymmetric" scalar channels (which are the co- 
sine and sine, respectively of the Wilson loop in question) 
and tensor glueball states. 

The wave function and mass are obtained from the 
correlation function: 



<?(»*;*) = <$t(f,t)*(o,o)), 



(2) 



where one needs to subtract the vacuum contribution 
from the correlator for ++ state. The source can be 
held fixed while the sink takes on the r dependence. This 
proves to be helpful in maintaining a good signal. The 
disentangling of the glueball and torelon is usually taken 
care of automatically by the choice of Wilson or Polyakov 
loops. 

To increase the overlap with the lowest state and re- 
duce the contamination from higher states, we exploit 
the APE link smearing techniques 0]. The procedure 
is implemented by an iterative replacement of the origi- 
nal spatial link variable by a smeared link. This results 
in correlations which reach their asymptotic behaviour at 
small time separations. In addition, the noise from ultra- 
violet fluctuations is reduced. The smearing parameter 
is fixed to 0.7 and ten iterations of the smearing process 
are used. To find the optimum smearing value, n, we 
examine the ratio (at r = and 1) 

C{r,t + 1)/C(r,t), 

which should be maximum for good ground state domi- 
nance. Using lxl loop as template, the best signal is 
obtained with four smearing steps, with lxl and 2x2 
loops being almost indistinguishable. At (3 = 2.0, the 



The link-link operator used in Ref. sums up a large number 
of loops; some of these loops have a zero winding number and 
project on glueballs - others have non-zero winding number and 
project on flux states also called torelons. 



signal in 1 x 1 showed a slow convergence with n, hence 
2x2 loops were preferred for optimum overlap. A typi- 
cal value which proved to be sufficient for this case was 
n = 2. 

A second pass was made to measure the optimized cor- 
relation matrices 

C«(t) = (3(r is t)*(r is 0)) - ^(n))^)). (3) 

Let ip^ be the radial wave function of the A:-th eigen- 
state of the transfer matrix, then 



(4) 



The glueball masses and the wave functions are extracted 
from the Monte Carlo average of CV,- (t) by diagonalizing 
the correlation matrices C(t) for successive times t: 



C{t) = R(t)D(t)R(t), 



(5) 



where D is a diagonal matrix of the eigenvalues and R 
a rotation matrix whose columns are the eigenvectors of 
C. Each eigenvector of C matches an eigenstate if)( k >(r) 
of the complete transfer matrix. As the wave function 
is largest at the origin, one would first determine the 
glueball mass with the optimal separation, and then fix 
that mass for all r, and extract the wave function for less 
optimal separations. Similar to the case of mesons |15| . 
the wave function is expected to decrease exponentially 
with the r at large separations and is therefore fitted with 
the simple form 



(6) 



to determine the effective radius Tq. The effective mass 
can be read off directly from the largest eigenvalue cor- 
responding to the lowest energy 



l eff 



= log 



Ao(r = 0,t=l) 



A (r = 0,t = 2) 



(7) 



III. SIMULATION RESULTS AND DISCUSSION 

A. Results at zero temperature 

Most of our Monte Carlo calculations are carried out 
on 16 2 x 16 lattice with periodic boundary conditions (16 2 
is the space-like box and 16 is the extension in Euclidean 
time direction). The gauge configurations are generated 
using the Metropolis algorithm. After the equilibration, 
configurations are stored every 250 sweeps; 2000 stored 
configurations are used in the measurement of glueball 
masses. Measurements made on the stored configura- 
tions are binned into 10 blocks with each block contain- 
ing an average of 200 measurements. The mean and the 
final errors are obtained using single-elimination jack- 
knife method with each bin regarded as an independent 
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data point. Three sets of measurements were taken at 
(3 = 2.0, 2.25 and 2.5. Some finite-size consistency checks 
are done at f3 = 2.25 2.5 on an 20 2 x 20 lattice. 

The glueball correlation function for the ++ channel 
against t at /3 = 2.25 is shown in Fig. ^ It can be 
seen that the expected behaviour of correlation function 
is attained virtually straight away. The absolute errors 
in the correlation functions are expected to be indepen- 
dent of t for large t. Our errors are consistent with this 
expectation. 



t = 1 and 2. We found a little variation (less than two 
percent) in the eigenvectors of C(t) with t which suggests 
that there is no mixing with states of distinct masses. 
Typical plots of the wave functions, normalised to unity 
at the origin, for the symmetric and antisymmetric scalar 
glueballs, at (3 = 2.0, 2.25 and 2.5 are shown in Figs. |3| 
and ^respectively. For guiding the eyes the Monte Carlo 
points of the same /3-value are connected with straight 
lines. The scalar wave function shows the expected be- 
haviour for all the (3 values analysed here. As for the 
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FIG. 1: Correlation function for the channel against t. 

Effective mass plot for f3 = 2.5 simulation is presented 
in Fig. For ++ and channels each it was possi- 
ble to find a fit region i mm — imax hi which convincing 
plateaus were observed. The effective masses are found 
stable using different values of t in Eq. 0, which sug- 
gests that the glueball ground state is correctly projected. 
At (3 — 2.5, we noticed considerable fluctuations in the 
tensor mass at large t. An acceptable fit was only possi- 
ble for t = [2-5]. 
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FIG. 2: Effective mass plot for scalar and tensor glueball 
states for j3 = 2.5. The dashed horizontal lines indicate the 
plateau values. 



FIG. 3: Scalar 0++ glueball wave functions measured on a 
16 3 lattice for various values of j3. 

antisymmetric channel we notice the presence of nega- 
tive contributions in the glueball wave function for r > 6 
at (3 = 2.5. However these contributions do not persist 
when the lattice size is increased from L = 16 to 20. This 
would mean that these effects are unphysical and can be 
described as a finite volume artifact. For this reason we 
extract the effective radius of the antisymmetric state 
from the wave function obtained at larger volume. The 
symmetric scalar glueball wave function, on the other 
hand, barely changes sign. 

Fig. shows the wave function for the tensor glueball, 
at (3 = 2.0, 2.25 and 2.5. The tensor wave function re- 
mains positive and shows the expected flatness. It can 
be seen that tensor glueball is much more extended than 
the scalar as one moves towards higher (3 values. This 
would imply that tensor is therefore more sensitive to 
the finite-size effects, which is very visible in the distor- 
tion of the wave function for large r at (3 = 2.5. Naively 
we would expect that the spatial size at which we begin 
to encounter large finite size effects to be related to the 
size of the glueball. 

The expected finite-size scaling behaviour of the mass 
gap near the continuum critical point in this model is 
not known; but Weigel and Janke [16J have performed 
a Monte Carlo simulation for an 0(2) spin model in 
three dimensions which should lie in the same univer- 
sality class, obtaining 



The Wave functions are extracted at time-separations 
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FIG. 4: Scalar glueball wave function on a 16 lattice 
at (3 = 2.0, 2.25 and 2.5. 



FIG. 6: The logrrithmic plot of the scalar and tensor glueball 
wave functions at /3 = 2.25 and 2.5. The effective radii can 
be obtained from the inverse slopes of the curves. 




FIG. 5: Tensor glueball wave function on a 16 lattice at 
(3 = 2.0,2.25 and 2.5. 



mates obtained with a x 2 /Ndf of 0.92 - 0.67. Due to dis- 
tortion 3 in the tensor wave function at small r at /3 = 2.5, 
a meaningful fit was possible only in the range 6 < r < 8. 
The effective radius obtained was confirmed by examin- 
ing the plateau in the ratio \og[\jj(r)/4>(r + 1)]. Note, 
that our lograthmically plotted wave functions (Fig. 0) 
are merely illustrations. 

To summarize: in the weak coupling region a spectrum 
of massive ++ , and 2 ++ glueballs is indicated with 

m(0—) < m(0++) < to(2++). 

Since there is a good signal for wave function persisting 
long enough to demonstrate convergence to the asymp- 
totic value, it seems to be reasonable to estimate mass 
ratios with our present method. The estimates of masses 
and ro, in lattice units, at various (3 values are shown in 
Tables HI IITI and ITTT1 



for the magnetic gap. In order to ascertain the finite-size 
effect on our measurements, we performed extra simula- 
tions on a 20 2 x 20 lattice at (3 — 2.25 and 2.5. The mass 
and size of ++ channel are almost unchanged as the lat- 
tice size increases from 16 to 20. We also find that our 
estimates for the tensor state are consistent with no fi- 
nite volume dependence at (3 = 2.25. However, the tensor 
mass was found to increase by about 4% and the effective 
radius by about 7% from 16 to 20 lattices at (3 = 2.5. We 
do not have enough data extrapolate mass and the radius 
to the infinite volume limit or to check whether the dif- 
ference is due to statistical errors or whether there is an 
incomplete convergence. Given that no mass reductions 
of sufficient magnitudes were found as the lattice volume 
is changed, none of our states could be interpreted as a 
torelon pair. 

In order to get some quantitative information on the ef- 
fective radius, the glueball wave functions are fitted in the 
range 3 < r < 8 by the form © . This form fits the data 
rather well for the scalar glueball with the best-fit esti- 



TABLE I: Masses of scalar glueballs in lattice units for two 
spatial extensions, L = 16 and 20. 







Mass 






0++ 




0" 




f3/L 


16 


20 


16 


20 


2.0 


0.803(6) 




0.441(4) 




2.25 


0.523(3) 


0.527(3) 


0.266(3) 


0.261(4) 


2.5 


0.364(3) 


0.369(2) 




0.182(2) 



Our results for lattice masses and mass ratios are gen- 
erally, within statistical errors, in agreement with the ex- 
isting Euclidean estimates 0, H3, LLa ! 11 perhaps a little 



3 Because of the distortion and impossible complete elimination of 
all the excited-states, especially near r ~ 0, it follows that Eq. 
IbI holds only in the limited interval, which does not include the 
vicinity of r ~ 0. 



TABLE II: Sizes of scalar glueballs in lattice units for two 
spatial extensions, L — 16 and 20. 





Size 




0++ 




-- 


/3/L 16 20 


16 


20 


2.0 1.39(16) 


1.01(12) 




2.25 2.75(43) 2.79(34) 


2.14(36) 


2.16(22) 


2.5 5.13(1.07) 5.20(94) 




5.05(1.05) 


TABLE III: Mass and size of tensor glueballs 


in lattice units 


for two spatial extensions, L — 


16 and 20. 




Mass 


Size 


/3/L 16 20 


16 


20 


2.0 1.18(1) 


5.02(68) 




2.25 0.82(2) 0.814(6) 


9.71(1.67) 


9.79(1.44) 


2.5 0.544(2) 0.583(1) 


10.13(2.56) 


10.19(2.37) 



high in places. Qualitatively our results, at zero temper- 
ature, are in agreement with the scenario of spectrum of 
massive magnetic monopoles [Tflj . 



TABLE IV: Glueball sizes in the units of string tension. 





K{= a 2 a) 






r —y/a 


r 2 ++y/o : 


2.0 


0.0508(5) 


0.0856 


0.31(14) 


0.24(9) 


1.13(18) 


2.25 


0.0221(3) 


0.0481 


0.40(17) 


0.32(16) 


1.14(22) 


2.5 


0.0119 


0.0272 


0.56(21) 


0.50(19) 


1.11(25) 



To extrapolate our effective radii to the continuum 
limit, we take the dimensionless products of sizes so that 
the scale, a, in which they are expressed cancels. We 
choose to take products of the effective radii, ro^/a, to 
a^/a since the string tension is our most accurately cal- 
culated quantity. As in the (3+l)D confining theories, 
we expect that dimensionless product of physical quanti- 
ties, such as ro^y/v, will approach their continuum limit 
with correction of 0(a%ff), where a e // is the effective lat- 
tice spacing in "physical units" when the mass gap has 
been renormalized to a constant 01 ■ The string tension, 
K(= a 2 a), is obtained by using the Wilson loop aver- 
ages and fitting the on-axis data with V(r). In Fig. [7| 
we show the product r^^y/a plotted against a e //- Since 
the products are plotted against a e f f , the continuum ex- 
trapolations are simple straight lines. We notice that the 
product ro z 2y/v varies only slightly over the fitting range. 
The non-zero lattice spacing values of the product are 
within 0.04 - 0.29 and 0.01 - 0.02 standard deviations of 
the extrapolated zero lattice spacing results for the scalar 
and tensor glueballs respectively. The striking feature of 
this plot is the little variation of the product with a e f t . 
This will make for very accurate and reliable continuum 
extrapolations. Linear extrapolations to the continuum 



limit yield values of 0.60 ± 0.05 and 1.12 ± 0.03, in the 
units of string tension, for the scalar and tensor states, 
respectively. In contrast to the tensor, the scalar glueball 
size shows significant finite-spacing errors. By setting the 
string tension to 420 MeV, we obtain the physical radii 
of 0.28(7) and 0.52(5) fm, for the scalar and tensor glue- 
balls, respectively. Our results show the size of the tensor 
glueball roughly two times as large as the scalar glueball. 
The extension of the method to the SU(3) lattice gauge 
theory yielded very similar results for the glueball 
size in the limit a — > 0. These estimates agree with the 
rough estimates of glueball sizes obtained at various tem- 
peratures in Ref. |9j. This is an improvement over the 
estimates obtained in Ref. Q where the predicted radius 
for the tensor glueball (~ 0.8 fm) was found four times 
larger than scalar glueball radius. 




FIG. 7: Glueball radii in the units of string tension as a 
function of the effective spacing, a e f / . Extrapolations to the 
continuum limit are shown as dashed lines. 



B. Finite temperature results 

To check the consistency of our method, we performed 
a study on an asymmetric lattice: 16 2 x 4 at (3 = 2.25. 
The procedure itself is a straight forward extension of the 
procedure adopted in the previous subsection. We do not 
plan to study the high temperature aspects of this model 
here but focus on the behaviour of the glueball mass and 
wave function in the deconfined region. 

The physical temperature T = l/(aN t ), is given via 
the lattice parameters as follows: 



1 



NtVK 



(9) 



For completeness, we give a temperature estimate of 
1.125 in the units of string tension. By setting the string 
tension to 420 MeV, we estimate a physical temperature 
of T ~ 1.25T C , where the T c ~ 360 MeV at pseudo- 
critical coupling fa = 1.87(2) 20]. One expects 



6 



that the high temperature phase has a massless photon 
and the linear potential is replaced by the two dimen- 
sional logarithmic Coulomb potential. This logarithmic 
behaviour is equivalent to a power-law dependence of the 
Wilson loop correlation function, 



C(r) = (Pt(r)P(O)) 



-V(T) 



(10) 



with an exponent which decreases as T increases. Fur- 
thermore, since the high-temperature phase of the gauge 
theory corresponds to the ordered phase of the spin sys- 
tem, the predicted power-law behaviour of the correla- 
tion function is just like that of a two-dimensional U(l)- 
invariant spin system - a 2-D XY model. 

Fig. |S| shows a plot of correlation functions versus sep- 
aration. The straight line indicates the fit to the form 
(|10|l . The finite temperature phase transition is visible 




log(r) 



FIG. 8: The logarithmic plot of the correlation function at 
[3 — 2.25. The straight line indicates power-law behaviour. 

in the change of the correlation function from exponen- 
tial to power-law behaviour. Thus it becomes evident 
that T > T c in our simulation. It can be seen that form 
(|10|l fits the data rather well. Nonetheless, our Monte 
Carlo simulations were unable to confirm that the expo- 
nent is mo ving towards the value of 0.25 (that of the 2-D 
XY model [23]) predicted for the continuum theory. Our 
estimated value for the exponent is four times larger than 
the predicted value. This indicates that our (3 value of 
2.25 is not large enough to gives us reason to hope that 
we are approaching continuum physics. An interesting 
feature to explore in this context is whether the cou- 
pling to the matter fields in the leading order ((3 — > 00) 
calculations will move the critical exponent towards the 
predicted value. 

Fig. shows the scalar and tensor wave functions ob- 
tained through the same analysis as in Figs. |21 and |S] 
Our results indicate that no significant changes occur in 
the scalar and tensor wave functions. Glueball masses 
appear almost unaffected. By comparing the results at 
T = and T = 360, we observe an effective mass reduc- 
tion, (am G (T - 0) - am G (T ~ 360)), of about 4%, with 
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FIG. 9: Scalar and tensor glueball wave functions measured 
on 16 2 x 4 lattice at /3 = 2.25. 



statistical uncertainties typically on less than a percent 
level, for ++ and 2 ++ glueball modes. This appear to 
be a very small change since we expect a rather continu- 
ous mass reduction of glueballs in the deconfined phase. 
This might be due to the fact that for zero momentum 
the power-law behaviour of the correlation function leads 
at short distances to the spin-wave results, which pre- 
vents us from seeing the massless excitations. The non- 
vanishing effective masses would suggest the presence of 
glueball modes above T c . Other work on finite tempera- 
ture su(3) um 

has also confirmed the survival of cor- 
relations above T c in the scalar and tensor colour-singlet 
modes. However, these studies have shown that ther- 
mal mass changes rather continuously across the critical 
temperature. The existence of the effective mass gives 
rise to the possibility that some of the nonperturbative 
effects survive in the deconfined phase, and the colour- 
singlet modes exist as metastable states above T c . The 
metastable states in the ordered phase (large (3) appear 
to be caused by the unusually large separation of a vortex 
pair, which may take many sweeps to recombine. Near 
the transition the number of vortices increase, and some 
of them begin to unbind. This eventually drives the sys- 
tem into a disordered phase as one moves to the region 
T <T C . Whether bound or metastable modes, the glue- 
balls can decay into two or more gluons thus acquiring 
finite width which is expected to become less negligible 
in the deconfined phase. Thus it becomes important to 
take into account the effect of width in best-fit analysis. 
This might also explain a very modest reduction of our 
masses at T > T c . However, from this study, it is not 
possible to determine whether such colour singlet modes 
really survive above T c as metastable states. An exten- 
sive systematic analysis, of unquenched improved lattice 
QCD at finite temperature, along these lines is under way 
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IV. SUMMARY AND CONCLUSION 

We have studied wave functions and sizes of scalar and 
tensor glueballs using improved 3-dimensional U(l) lat- 
tice model. In this preliminary study we take a more 
direct approach to the problem; instead of fixing a gauge 
or a path for the gluons, we measure the correlation 
functions from our lattice operators from spatially con- 
nected Wilson loops which, being the expectation values 
of closed- loop paths, are gauge invariant. This approach 
has the advantage that the disentangling of the glueball 
and torelon is usually taken care of automatically by the 
choice of Wilson or Polyakov loops. We observed that 
the size of tensor glueball is roughly two times larger 
than the size of the scalar glueball. We have successfully 
extended our method to extract the results of glueball 
sizes for SU(3) lattice gauge model and obtained some 
encouraging results |19f. We believe that our estimates 
are more reliable than the results obtained in Ref. Q, 
where the size of the tensor glueball was found to be ~ 0.8 
fm, four times as large as its scalar counterpart. The pre- 
dicted zero lattice spacing results are not actually found 
by extrapolation to zero lattice spacing, but are obtained 
instead from calculations at (3 of 2.2 of glueball size, with 
no accurate representation of the effect of the absence of 
extrapolation. Also the results were of limited interest 



because of their manifest dependence on the gauge cho- 
sen and the problem of disentangling of the glueballs and 
torelons. 

Finally, for completeness, we extended our method to 
measure the wave function and mass for a finite temper- 
ature deconfinement phase. For the lowest ++ and 2 ++ 
glueballs, no significant mass reduction was observed in 
deconfincd phase, while the wave functions remain almost 
unchanged. The existence of the effective mass indicates 
that colour-singlet modes may survive in the deconfincd 
phase as metastable states. In such a case glueball de- 
cay and decay width, as spectral component, in the de- 
confinement phase are the most feasible candidates for a 
more reliable analysis for the future studies. 
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